Nonlinear theory of resonant slow waves in anisotropic and dispersive plasmas 
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The solar corona is a typical example of a plasma with strongly anisotropic transport processes. 
The main dissipative mechanisms in the solar corona acting on slow magnetoacoustic waves are 
the anisotropic thermal conductivity and viscosity. Ballai et al. [Phys. Plasmas 5, 252 (1998)] 
developed the nonlinear theory of driven slow resonant waves in such a regime. In the present 
paper the nonlinear behaviour of driven magnetohydrodynamic waves in the slow dissipative layer 
in plasmas with strongly anisotropic viscosity and thermal conductivity is expanded by considering 
dispersive effects due to Hall currents. The nonlinear governing equation describing the dynamics 
of nonlinear resonant slow waves is supplemented by a term which describes nonlinear dispersion 
and is of the same order of magnitude as nonlinearity and dissipation. The connection formulae are 
found to be similar to their non-dispersive counterparts. 

PACS numbers: 52.20.-j; 52.25.Fi; 52.30.Cv; 52.35.-g; 96.60.P- 



I. INTRODUCTION 

Resonances are ubiquitous everytime magnetohydro- 
dynamic (MHD) waves are driven in a transversally in- 
homogeneous (relative to the ambient magnetic field) 
plasma. From a mathematical point of view a resonance 
is equivalent to regular singular points in the equations 
describing the dynamics of waves, but these singularities 
can be removed by, for example, dissipation. At res- 
onance interacting dynamical systems can transfer en- 
ergy to each other. In this context, resonant absorption 
has been suggested as a method to create supplementary 
heating of fusion plasmas, but was, however, later re- 
jected due to technical difficulties (see, e.g. Refs. [H-[3|)- 
lonson suggests, for the first time, that resonant MHD 
waves may be a means to heat magnetic loops in the solar 
corona. Since then, resonant absorption of MHD waves 
has become a popular and successful mechanism for pro- 
viding some of the heating of the solar corona (see, e.g. 
Refs. @, 01' ^^'^ references therein). More recently res- 
onant absorption has acquired a new applicability when 
the observed damping of waves and oscillations in coro- 
nal loops has been attributed to resonant absorption. 
Hence, resonant absorption has become a fundamental 
constituent block of one of the newest branches of solar 
physics, called coronal seismology (see, e.g. Refs. 

M)- 

A driven problem for resonant MHD waves occurs 
when there is an external (or internal) source of energy 
that excites the plasma oscillations. If there is a small 
amount of dissipation present in the system, after some 
time the system will attain a steady state in which all per- 
turbed quantities will oscillate with the same frequency 
uj. In a weakly dissipative plasma, dissipation only be- 
comes important in thin layers enclosing the singularity 
(the so-called dissipative layer). Two types of driving 
are possible. First, direct driving, where the source of 
energy is inside the system or at the boundary and the 
oscillations of resonant magnetic field lines are directly 
excited by this source of energy. In the context of res- 



onant absorption, this type of driving was studied by, 
for example, Ruderman et al. [ll|; and Tirry et al. 
[T^ ]. Secondly, in the case of indirect or lateral driv- 
ing, the energy source can be either outside or inside the 
system. This energy source excites fast or slow magne- 
tosonic waves which propagate across and along magnetic 
surfaces and reach the resonant magnetic surface where 
their energy is partly dissipated due to resonant coupling 
with localized Alfven or slow waves. The lateral driving 
problem was studied by, for example, Davila for pla- 
nar geometry and by Erdelyi pisj for cylindrical geome- 
try. In the present paper, we consider only the lateral 
driven case (for a comprehensive background to lateral 
driving see, e.g. Refs. [l3|-|2^)- An important property 
of these waves is that their damping rate is almost in- 
dependent of the values of the dissipative coefficients, a 
situation characteristic for dissipative systems with large 
Reynolds numbers. As a result, the damping rate of near 
resonant MHD waves can be many orders of magnitudes 
larger than the damping rate of MHD waves with the 
same frequencies in homogeneous plasmas. The damp- 
ing allows the waves energy to be converted into heat, 
which has made resonant absorption a subject of intense 
study. 

The concept of connection formulae was introduced by 
Sakurai et al. (2]| . to determine the jumps in the normal 
component of velocity and perturbation of total pressure 
across the dissipative layer. The connection formulae 
avoid solving the full dissipative MHD equations when 
studying resonant MHD waves. Instead, the narrow layer 
embracing the resonant point can be considered as a sur- 
face of discontinuity and ideal MHD equations can be 
used to the left and right of this surface. Connection for- 
mulae are utilized in order to connect the ideal solution 
over the discontinuity. 

Most studies on driven resonant MHD waves use 
isotropic viscosity and/or electrical resistivity. However, 
the solar corona is a well-known exarnple of a plasma 
where viscosity is strongly anisotropic |22l | . HoUweg and 
Yang [2^ studied the laterally driven problem in the cold 
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plasma approximation. They found that anisotropic vis- 
cosity does not remove the Alfven singularity (if the Bra- 
ginskii's viscosity tensor is approximated by its first term 
only). If Braginskii's full viscosity tensor is considered, 
the Alfven singularity would be removed by the shear vis- 
cosity. The way the dissipative term appears in the gov- 
erning equation is, however, identical to that of isotropic 
viscosity. 

For the case of slow resonant waves the situation is 
different. The laterally driven linear slow resonant waves 
in plasmas with strongly anisotropic viscosity and ther- 
mal conductivity was studied first by Ruderman and 
Goossens '2^ . They successfully showed that anisotropic 
viscosity and/or thermal conductivity removes the slow 
resonance present in ideal plasmas. They also obtained 
the explicit connection formulae, which are identical to 
those found in the case of plasmas with isotropic viscosity 
and finite electrical resistivity, [2l|. This fact supports 
the hypothesis that in weakly dissipative plasmas the 
connection formulae are independent of the exact form 
of dissipative processes present in the dissipative layer. 

The laterally driven nonlinear slow resonant waves in 
plasmas with strongly anisotropic viscosity and thermal 
conductivity was first studied by Ballai et al. [1^ . They 
found that nonlinearity was crucial in the dissipative 
layer. The governing equation for slow wave dynamics 
in the dissipative layer was derived and the implicit con- 
nection formulae were found, (explicit connection formu- 
lae have only been found for linear theory and for the 
limit of strong nonlinearity, [1^). The governing equa- 
tion was almost identical to that found by Ruderman et 
al. [13], however the dissipative term was laterally de- 
pendent (9) rather than normally dependent (^) depen- 
dent. The implicit connection formulae found coincide 
with those found in plasmas with isotropic viscosity and 
finite electrical resistivity [l^] ■ 

A drawback of previous studies on resonant absorption 
is that even though anisotropy is considered, dispersion 
(by, e.g. Hall effect) is neglected. This approximation 
is acceptable only for lowest regions of the solar atmo- 
sphere. The solar corona is known to be strongly mag- 
netized, hence the Hall term can be comparable with 
other effects considered in the process of resonance. The 
present paper will extend the nonlinear theory of res- 
onant slow MHD waves in the dissipative layer with 
strongly anisotropic viscosity and thermal conductivity 
to include Hall dispersion and show that the effect of 
this new addition is of the same order of magnitude as 
nonlinearity and dissipation near resonance. The paper 
is organised as follows. In the next section we introduce 
the fundamental equations and discuss the main assump- 
tions. In section HI we derive the nonlinear governing 
equation for waves in the dissipative layer. Section IV is 
devoted to the derivation of the nonlinear analogues of 
the connection formulae. Finally, in section V we summa- 
rize and draw our conclusions pointing out a few further 
applications to be carried out in the future. 



II. FUNDAMENTAL EQUATIONS 

In what follows we use the visco-thernial MHD equa- 
tions with strongly anisotropic viscosity and thermal con- 
ductivity. We assume that the plasma is strongly magne- 
tised, so that the conditions oJeTe » 1 and uiin ^ 1 are 
satisfied, where 0Je{i) is the electron (ion) gyrofrequency 
and Te(j) is the mean electron (ion) collision time. Due 
to the strong magnetic field, transport processes are de- 
rived from Braginskii's expression for the viscosity tensor 
TT (see Refs. and [1^). Under coronal conditions it 
is a good approximation to retain only the first term of 
Braginskii's expression for viscosity, [22| . namely 

n^?Jo(h^h- i/^ [3b • (b • V)v - V • v] , (1) 

where v is the velocity, h — H/B is the unit vector in the 
direction of magnetic field and rfo is the first Braginskii 
coefficient of viscosity (compressional viscosity) . / is the 
unit tensor and the symbol (g) denotes the dyadic product 
of vectors. 

In a strongly magnetised plasma the thermal conduc- 
tivity parallel to the magnetic field lines dwarfs the per- 
pendicular component so the heat flux becomes, [2§|, 



-K||b(b- VT), 



(2) 



where T is the temperature and K|| is the parallel co- 
efficient of thermal conductivity (in the solar corona, 
Kji = 9 X 10-122^5/2 Wm-^K-i). 

In the solar corona the finite electrical resistivity can 
be neglected as it is several orders of magnitude smaller 
than the dissipative coefficients considered here, see e.g. 

Using Eqs. ([T]) and ([2]), the visco-thermal MHD equa- 
tions are 



0, 



(3) 



- + (v.V)v 



-ivP+ — (B • V)B 
P PP 

+ i(V-b) X (,7obg), (4) 
P 



-.Vx(vxB) 



dT — — 

— V • VT + (7 - i)rv • V = 
at 



H, 



(5) 



7-1 
Rp 



|v- [S||b(b.VT)] +i^oQ'|, (6) 



5- + ^Q, 
2fi 3 ^' 



(7) 
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P = RpT, 



V • B = 0, 



g = 3b • (b • V)v - V ■ V. 



(8) 



(9) 



(10) 



In Eqs. ([5|)- pU|) . p is the kinematic pressure, p the den- 
sity, P the viscosity modified total pressure (kinetic and 
magnetic), 7 the adiabatic exponent, R the gas constant 
and p, the magnetic permeabihty. The term H in Eq. ([5]) 
is the Hall term given by 



H = —V X f — B X V X B 



(11) 



where e is the electron charge and rig is the electron num- 
ber density. The propagation of compressional linear and 
nonlinear MHD waves in Hall plasmas has been studied 
by, for example, Baronov & Ruderman ^31j, Ruderman 
[Si; [13; m, Ballai et al. ^ and Miteva et al. [H]. As 
stated in an earlier study by Huba [ST*! , Hall MHD is only 
relevant to plasma dynamics occurring on length scales 
of the order of the ion inertial length, di — c/uoi^ where c 
is the speed of light and uji is the ion plasma frequency. 
For the present paper this would require that di ~ ^{Sc), 
where Sc is the thickness of the dissipative layer. Indeed, 
starting from the upper chromosphere this condition is 
satisfied and the lengths involved in the problem are of 
the order of 10 — 100m. Hall currents arise when consid- 
ering off diagonal terms in the conductivity tensor. 

We adopt a Cartesian coordinate system, and limit our 
analysis to a static background equilibrium (vq = 0). We 
assume that all equilibrium quantities depend on x only. 
The equilibrium magnetic field, Bq, is unidirectional and 
lies in the yz-plane. The equilibrium quantities must 
satisfy the condition of total pressure balance. 



^0 

Po H = constant. 

2fi 



(12) 



In addition we assume that the wave propagation is inde- 
pendent ofy {d/dy = 0). In linear theory of driven waves 
all perturbed quantities oscillate with the same frequency 
Lu, SO they can be Fourier-analysed and taken to be pro- 
portional to exp(i[kz — ujt]), so solutions are sought in the 
form of propagating waves. All perturbations in these 
solutions depend on the combination 9 = z — Vt, rather 
than z and t separately, with V = uj/k. In the context of 
resonant absorption the phase velocity, V, must match 
the projection of the cusp velocity, ct, onto the z-axis 
when X = Xc- To define the resonant position mathemat- 
ically it is convenient to introduce the angle, a, between 
the z-axis and the direction of the equilibrium magnetic 
field, so that the components of the equilibrium magnetic 
field are 



The definition of the resonant position can now be writ- 
ten mathematically as 



V = Ct (xc) cos a. 
The square of the cusp speed is defined as 



Ct 



(14) 



(15) 



where the squares of the sound and Alfven speeds are 
given by 



IPo 



Po MPo 



(16) 



with all speeds being dependent on the coordinate x. 

In a nonlinear regime the perturbations cannot 
Fourier-analysed, however, in an attempt to adhere as 
close to linear theory as possible we look for travelling 
wave solutions and assume all perturbed quantities de- 
pend on 9 ^ z — Vt where V is given by Eq. (fT4|) . The 
perturbations of the physical quantities are defined by 

p = po + p, p = po+p, T = To+T 

B = Bo + B, H = Ho + H, P = Pq + P (17) 

It is clear from Equations (jlip and (fT7| that Hq = 0. 

The dominant dynamics of resonant slow waves, in lin- 
ear MHD, resides in the components of magnetic field 
and velocity that are parallel to the equilibrium mag- 
netic field (as well as in the compressional quantities p, 
p and T). This dominant behaviour is created by an 
x~^ singularity in the spatial solution of these quantities 
at the cusp resonance, (jlH); these are known as large 
variables. The normal components of velocity, u, and 
magnetic field, B^, are also singular however their singu- 
larity is proportional to In |a;|. In addition, the quantities 
P and the components of B and v that are perpendicular 
to the equilibrium magnetic field are regular; these are 
known as small variables. 

To make the mathematical analysis more concise and 
to make the physics more transparent we define the com- 
ponents of velocity and magnetic field that are parallel 
and perpendicular to the equilibrium magnetic field (as 
well as existing in the yz-plane): 



v± 



I - 



By Bz 



V — w 

By - Bz 



sma 
cosa 



cos a 
sin a 



(18) 



Oy 



i?osina, B{ 



Oz 



Bq cos a. 



(13) 



where u, w, By and Bz are the y- and z-components of 
the velocity and perturbation of magnetic field, respec- 
tively. 

Let us introduce the characteristic scale of inhomo- 
geneity, Imh- The classic Reynolds number, i?e, and the 
Pecklet number, P^, are defined as 

_ Pa^Vknh p _ PO^RVUnh ^-^^^ 
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These two numbers determine the importance of vis- 
cosity and thermal conduction. We introduce the total 
Reynolds number as 



1 
R 



1 

Re 



1 



(20) 



The aim of this paper is to derive the governing equa- 
tion for waves in the slow dissipative layers taking into 
account nonlinearity, dissipation and dispersion simulta- 
neously. Since for the coronal plasma i? ^ 1, we con- 
sider only weakly dissipative plasmas. We introduce, e, 
as the dimensionless amplitude of variables far away from 
the resonance. Linear theory predicts that the charac- 
teristic scale of dissipation, Idiss, is of the order R~^linh- 
The typical largest nonlinear term in the system of MHD 
equations is of the form gdg/dz while the typical dissipa- 
tive term if of the form fj^d^g/dz^, where g is any 'large' 
variable. Linear theory shows that 'large' variables have 
an ideal singularity {x — Xc)~^ in the vicinity of x = Xc- 
This implies that the 'large' variables have dimensionless 
amplitudes, inside the dissipative layer, of the order of 
eR. It is now straightforward to estimate the ratio of a 
typical nonlinear and dissipative term. 



gdg/ dz 
f]od'^g/dz'- 



(eR:' 



(21) 



where the quantity (j) can be considered as a nonlinearity 
parameter. If the condition eR^ <C 1 is satisfied, linear 
theory is applicable. On the other hand, if eR^ > 1 then 
nonlinearity has to be taken into account when study- 
ing resonant waves in the dissipative layers. Using the 
same scalings, Ballai et al. 25] showed that nonlinearity 
has to be considered whenever slow wave resonant ab- 
sorption is studied in the solar corona. This conclusion 
simply means that in the solar upper atmosphere reso- 
nant absorption is a nonlinear phenomena. 

In order to obtain nonlinearity and dissipation of equal 
order we assume that eR^ — ^{1), i.e. R ^ e^^^^, when 
deriving the nonlinear governing equations for waves in 
the slow dissipative layer. Accordingly, we can re-scale 
the coefficients of viscosity and thermal conductivity as 



'70 



(22) 



du dP Bo cos a dB^ 



du 



de =PVd-x^'"m 



pV 



du dB^ B^ dB^ 



36 



fi dx p, 89 



^"^Yx^^ + ^Q^Ai^lM), (25) 



§9 I ^"^"^ 



P sin a 



Bo cos a 



P u 



dx 



dv± dv± \ 

I. 1- 7/7 1 



B^ dBi_ 



B, dBj 



fi dx fi 89 



9 / ~ Bo cos a 
poVv\\ — i^cosa-l B|| 



89 



dvn 



dv 



dv\\ 



Bx dBo 



/i dx 

Bx dB\\ B, dB\\ 



p dx p 89 



VBx + Boucosa — wB^ — uBz 



1/2 ^-^11 

e' x~:^ cos a sin a, (28) 
o9 



d ,,,r, o ^ d{uBx_) d{wB^) 



B 



ox 



" 89 ^8x89 



cos a, (29) 



89 



{VB\\ + BqV\\ cos a 



8{Bou) ^ 8w 

8x ° 89 

8v 



8{uB\\) ^ 8{wB^) ^ dv^ 
8x 89 ^ 8x 



B, 



II 

89 



1/2 X 9B\\ 8p . 



po 89 8x 



sin a, (30) 



We also consider the coefficient of Hall conduction, de- 
fined as X = TjujeTe, where rj is the coefficient of mag- 
netic diffusivity (although ry is small enough, in the solar 
corona, to be neglected in comparison to 770, here it is 
multiplied by the product UeTe which is very large under 
coronal conditions). Similar to the previous dissipative 
coefficients, we introduce the scaling 



X 



(23) 



Using Eq. (f22|) . we can rewrite Eqs. ([3])- ([7]) in the form 



U^ 
89 



8{pou) 
8x 



Po 



8w 
89 



8{pu) 8{pw) 



8x 



89 



(24) 



8Bx 8B, 



8x 



89 



0, 



,,8T dTo fdu 8w\ 

8T 8T , , f8u dw 



1/2 7-1 / I 2 



8x 



XKII 



89^' 
dT 



^ dx 8x J ^89 



(31) 



(32) 
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P 



B1 +Bi 



2BqB\ 



du 
dx 



Top - poT = pT, 



(33) 



(34) 



du 
'80 



36,1 h 



' dx 
du 
dx 



dv\\ 
'■~d0 



dw 



(35) 



We should state that in Eqs. we have used 

the coefficient of Hah conduction, x, which does not 
contribute to the total Reynolds number because it is 
the multiplier of dispersive terms rather than dissipative 
ones. The derivation of the expressions of the Hall terms 
in the induction equation can be found in the Appendix. 

The set of Eqs. (I111)-(I23, will be used in the next 
section to derive the governing equation for wave motion 
in the dissipative layer. 



III. DERIVING THE GOVERNING EQUATION 
IN THE DISSIPATIVE LAYER 

In order to derive the governing equation for wave mo- 
tions in the slow dissipative layer we employ the method 
of matched asymptotic expansions (J38|). This method 
requires us to find the so-called outer and inner expan- 
sions and then match them in the overlap regions. This 
nomenclature is ideal for our situation. The outer expan- 
sion corresponds to the solution outside the dissipative 
layer and the inner expansion corresponds to the solution 
inside the dissipative layer. A simplified version of the 
method of matched asymptotic expansions, developed by 
Ballai et al. is adopted here. 

We only consider weakly dissipative plasmas, so vis- 
cosity and thermal conductivity are only important in 
the narrow dissipative layer (here dissipation and non- 
linearity are of the same order) embracing the resonant 
position. Far away from the dissipative layer the ampli- 
tudes of perturbations are small, so we use linear ideal 
MHD equations in order to describe the wave motion far 
away from the dissipative layer. The full set of nonlinear 
dissipative MHD equations are used for describing wave 
motion inside the dissipative layer where the amplitudes 
are much larger than those far away from the dissipa- 
tive layer. We therefore look for solutions in the form 
of asymptotic expansions. The equilibrium quantities 
change only slightly across the dissipative layer so it is 
possible to approximate them by the first non- vanishing 
term in their Taylor series expansion with respect to x. 
Similar to linear theory, we assume the expansions of 
equilibrium quantities are valid in a region embracing 
the ideal resonant position, which is assumed to be much 



wider than the dissipative layer, hnh/ldiss R- This im- 
plies there are two overlap regions, one to the left and 
one to the right of the dissipative layer, where both the 
outer (the solution to the linear ideal MHD equations) 
and inner (the solution to the nonlinear dissipative MHD 
equations) solutions are valid. Therefore, both solutions 
must coincide in the overlap regions which provides the 
matching conditions. 

Before deriving the nonlinear governing equation we 
ought to make a note. In linear theory, perturbations of 
physical quantities are harmonic functions of 9 and their 
mean values vanish over a period. In nonlinear theory, 
however, the perturbations of variables can have non- 
zero values as a result of nonlinear interaction of differ- 
ent harmonics. Due to the nonlinear absorption of wave 
momentum, a mean shear flow is generated outside the 
dissipative layer, as shown by Ruderman et al. [27| in 
Cartesian geometry. In our scenario a mean shear flow 
is created outside the dissipative layer, but as there is no 
perpendicular component to viscosity oscillating plasma 
can slide past each other without friction. This produces 
a mean flow with infinite amplitude. However, bound- 
aries can prevent such oscillations. Therefore, we will 
assume such boundaries exist and will not consider the 
generation of mean shear flow. 

The first step in our mathematical description is the 
derivation of governing equations outside the dissipative 
layer where the dynamics is described by ideal (t^o = 
K|| = 0) and linear MHD. Since dispersion is assumed to 
act over scales comparable to dissipative scales, disper- 
sion will only be taken into account inside the dissipative 
layer. The linear form of Eqs. ((24l) - ((35)l can be obtained 
by assuming a regular expansion of variables of the form 



(36) 



and collect only terms proportional to the small param- 
eter e. This leads to a system of linear equations for the 
variables with superscript '1'. All variables can be elimi- 
nated, with the exception of u'-^-' and P^^\ leading to the 
system 



Po 



Adu(^ 



where 



dx 



F 



F de 



dx 



PoAC 



V de 



y4 _ y2 



cos^ a 



A^V^ -vl cos^ a. 



4) 



cos^ a) 



(37) 



(38) 



(39) 



The reason for using P instead of P is that outside the 
dissipative layer the plasma is ideal, so there is no viscous 
addition to the total pressure. The remaining variables 
can all be expressed in terms of u^^^ and P^^\ 



,(1) 



T^sina 
PoA 



„(i) 



Vcl cos a 
PoC ' 



(40) 



6 



(1) _ Bp cos a _ cos g sing 



PoA 



asj'^ Bo „ c| cos2 a) ap(i) dBo 



de 



PoC 



de V dx 



dx 



de 



de 



:s 

c de 

V^dP^ 

~ "c de 



,(1) 



V dx '■ 



dTW _ l)Tol/2 9p(i) 



961 



961 



(41) 



, (42) 



(43) 



(44) 



(45) 



-^VR dx 

The differential Eq. ([57)1 have regular singularities at 
the resonance, therefore the solutions can be obtained in 
terms of Frobenius series with respect to x (for details 
see, e.g. [Hi ; [131) of the form 

p(i) = pa)(0) ^ p!^^\e)x\n \x\ + p^^\e) + . . 



,(1) = 



u\'' [e] In |a;| + u^^'i^) + Uz'iS)^ In \x\ 



(1) 



, (46) 
.. (47) 

The coefficient functions depending on e in the above 
expansions are generally different for a; < and x > 
0. The salient property of these solutions is that the 
perturbation of the total pressure is regular at the ideal 
resonant position x — Xc- From Eqs. (|10)) - (|45p . we see 
that the quantities v'^^ and i?^^"* are also regular at a; = 
Xc- AH other quantities are singular. The quantities u^^^ 
and b'^x'^ behave as In \x\, while ^y^^' P'''^^ ^"^^ 
T'^^^ behave as x~^ , so they are the most singular. This 
property can be extended and will hold to all higher order 
approximations, [27j . 

Now let us concentrate on the solution in the dissi- 
pative layer. The thickness of the dissipative layer is of 
the order ImhR^^- We have assumed that R ~ e^^/^ 
so we obtain hnhR^^ — ^i^^^^linh)- The implication of 
this scaling is that we have to introduce a new stretched 
variable to replace the transversal coordinate in the dis- 
sipative layer, so we are going to use ^ = e~^^'^x. Equa- 
tions ((24l) - ((35l) are not rewritten here as they are easily 
obtained by the substitution of 



(48) 



A = £-1/2 A 

dx 9^' 

for all derivatives. The equilibrium quantities still de- 
pend on x, not ^ (their expression is valid in a wider 
region than the characteristic scale of dissipation). All 
equilibrium quantities are expanded around the ideal res- 
onant position, X = Xc, as 



/o = /o. + e 



dfo 
d^ 



dfo 
dx 



(49) 



where /o is any equilibrium quantity and the subscript 
'c' indicates the equilibrium quantity has been evaluated 
at the resonant point (we can always make = by 
proper translation of the coordinate system) . 

We seek the solution to the set of equations obtained 
from Eqs. (l24])-([35]) by the substitution of x = e^/^^ into 
variables in the form of power series of e. These equations 
contain powers of e^^^, so we use this quantity as an 
expansion parameter. To derive the form of the inner 
expansions of different quantities we have to analyze the 
outer solutions. First, since v± and B± are regular at a; = 
Xc we can write their inner expansions in the form of their 
outer expansions, namely Eq. p6p . The quantity P is 
the sum of the perturbation of total pressure P, which is 
regular at x = Xc, and the dissipative term proportional 
to Q. From Eq. ((35|) it is obvious that Q behaves as x~^ 
in the vicinity of a; = Xc- Far away from the dissipative 
layer, Q is of the order e. Since the thickness of the 
dissipative layer is of the order e^^'^knh, Q is of the order 
e^^'^ in the dissipative layer. However, Eq. ([33|) clearly 



shows the term proportional to Q contains a multiplier, 
e^/^, which implies the contribution of P supplied by the 
dissipative term is of the order e. As a consequence, we 
write the inner expansion of P in the form of its outer 
expansion, Eq. ((36)) . The amplitudes of large variables 
in the dissipative layer are of the order e^^"^, so the inner 
expansion of the variables t;|| , _B|| , p, p and T is 



£5 



(2) 



(50) 



The quantities u and B^ behave as ln|a;| in the vicin- 
ity of a; = a;c, which suggests that they have expansions 
with terms of the order of elne in the dissipative layer. 
Ruderman et al. [2^ showed that, strictly speaking, the 
inner expansions of all variables have to contain terms 
proportional to e In e and e"^/^ In e. In the simplified ver- 
sion of matched asymptotic expansions, [25l |. we utilize 
the fact that | lne| <C e~'^ for any positive k and e — > -1-0, 
and consider In e as a quantity of the order of unity. This 
enables us to write the inner expansions for u and B^ in 
the form of Eq. (|36| . 

We now substitute the expansion ([36]) for u, B^, P, 
v±, B± and the expansion given by ()50|) for v\\, _B||, p, p, 
T into the set of equations obtained from Eqs. (|24p - (P5)) 
after substitution of a; = e^/^^. The first order approxi- 
mation (terms proportional to e) , yields a linear homoge- 
neous system of equations for the terms with superscript 
'1'. The important result that follows from this set of 
equations is that 



p(i) =pW(0), 



(51) 



that is to say P*^^^ does not change across the dissipa- 
tive layer. This result parallels the result found in lin- 
ear theory (^; ^39]) and nonhnear theory ([111; [27j'). 
Subsequently, all remaining variables can be expressed in 
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terms of u^^\ v\P and P^^^ as 



,(1) _ ^sj^pn)^^^^ ^(1) ^ _ Bo, cos a 



Bq^ Cg^ sin a cos a 



u^^\ (52) 



(53) 

i?r = - = (54) 



(1) _ Bq^V ji) _ Pa.Vji) 



cos a 



(1) ^ _PoJ;^^(i) yd) ^ (7-1)^0,^ ^(1)^ 



2 I 
c^^ cos a ' 



(55) 



In addition, we find that the equation that relates u^^^ 



and v'^^^ is 



+ 



y2 



= 0. 



(56) 



In the second order approximation we use only the 
expressions obtained from Eqs. 

(pi)) and ([33]). Employing Eqs. ([51]) -((561), we replace the 
variables in the first order approximation. The equations 
obtained in the second order approximation are 



ay(2) auf) \ ap(2) 

POe I — I COS a I - V- 



80 



80 



V 



80 " Wa; 



dBo 
dx 



dv 



(1) 
I 



V{v 



39 



Bo,c|^ sin^ a dP^'^^ 



' " cos a T^v;, 



d0 

(i)~ 
I 



Bo^ sin a 



''Ac 



ci 



(59) 



,(1) 



961 



(7 " l)To. 



^1 +<^-') 



sin^ a dP 



961 



cos a 



(7 — l)y'«|| COS a a^U|| 



a (1) 
OTii 

56* 



i(i) az;f|'^ 



7^^A„ 



I COS a 8£, 



4. (1)^^ 

^t^M 



96*2 
(i)~ 



iZk„(2) 



(2) 



da; 



+4, 56* 
_^ (7 - 1)^0, / dpo 

(7 - i)ro,poe^i, 
4. (^^1 + 4J 



(60) 



(1) 



(«!'^)', (61) 



V 



0%^ sin^ a dP po.V ( u^^^ ^^^f 



(1) 



d0 



COS a 8S, 



(1) 



■^^^.^Wr^ I , (57) 



96* 



(2) , ^Oc 0(2) POc^ /rf5o\ . (1) 

II i?Oe COS a \ dx J " 



1 



34. + 2t'i, 



(1) 



poc4e 



-770 cos a — 5 5 — —^^5 5- 

3 + 96* 2 + 



("1")' 



(62) 



In deriving the above system we have used the fact that 



'Scl + 2.1 
^ - .2 +,2 -^'^osa, 



(63) 



Vv 



(2) 



-Boe COS a (2) 
i>ii 



MPOe 



i?o COS a 
H ^ 1 

Pl^POc 

POc V d-x 



XI) 



dBo 
dx 



cos a 

dpW 



POc 

V 



Vo 



d0 

dBo ' 
dx 

2„(1) 



cos2 a 2vl^ + 3c|^ d'v\ 



POc 



902 



(58) 



obtained from Eq. psp in the first order of approxima- 
tion. With the exception of Eq. ([59]) . which has the 
addition of the Hall term, these equations are identical 
to those found by Ballai et al. [25^ . 

The left-hand sides of the set of Eqs. could 
be obtained from the left-hand sides of the first order 
approximation by substituting variables with the super- 
script '2' for those with superscript '1'. The first order 
of approximation possesses a non-trivial solution, so Eqs. 
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(l57|) - (|62p are compatible only if the right-hand sides of 
Eqs. ([57)) - ((^ satisfy a compatibility condition. To de- 
rive the compatibility condition we express p*^^-' and b'"^'^ 

in terms of m'-^-', u^^\ v'^^^ and P^^\ using Eqs. ([55)1 
and ((5n|) - ([5^ . Subsequently, we substitute these expres- 
sions into Eqs. (|57p and ((59)) . to obtain 



cos a 89 



V{i 



c| sin^ a) 



V 



POc'^Ac ^^^^ ^ 



de 



V a cos a 



2 /dSo 



?7oVcosa f'2v\^+3cg^ 



1 



(1) c'^ll 
II 89 



,(1) 



de 

8vlP 



cos a 

(64) 



Xsina ^''^[i ' ^^\\ "* 



where we have used the notation 
_ ^ [(7 + 1)^1 + 3c|J ^2 ^ cosa 

Oj — 2 ) 



(67) 



7P0c^4, 



Si = ^ — cos asm a, 



dx 



(69) 
(70) 



Equation (I66p differs from its counterpart found by Ballai 
et al. [25] only by the last term of the left-hand side 
representing Hall dispersion. If we linearize Eq. ([66|) 
and take v''^^ proportional to exp(z/c0), we arrive at the 
linear equation for the parallel velocity obtained in linear 
theory by Ruderman and Goossens [2J]. 

Equation (j66p is the complete nonlinear governing 
equation for the parallel velocity in the slow dissipative 
layer. The function P^^^ in this equation is determined by 
the solution outside the dissipative layer, and is thought 
to be the driving term. We should note here that the 
dispersion (the last term on the left-hand side) appears 
as a nonlinear term {nonlinear dispersion). 



a (2) 

at; 1 1 



9^ v\ cosa 89 



sm a 



P0.^w2 



1/2 



To^c|^ cosa 



dx 



dv 



(1) 



d9 
V 



89 



dv\\ 



Cg cosa 9^ 



Ki24,±il±l)!H„(i)^ 



(1) 



l^cosa 

IPO A. 



27% /2t; 



3ci 



+ 



II 89 



R 



.2„(1) 



89^ 



(65) 



It can be seen that Eqs. ((64|) and ((65l) have identical left- 
hand sides. Extracting these two equations we derive the 
compatibility condition, which is the equation connecting 
v\^^ and P(i) 



AC- 



dv 



(1) 



,(i)^^l 



(1) 



1/3 A 



de 



89 



de^ 



n- 



Vcg^ COS a 



(66) 



IV. NONLINEAR CONNECTION FORMULAE 

In linear dissipative MHD it was assumed that when 
dissipative effects are weak they are only important in 
the thin dissipative layer that embraces the ideal reso- 
nant position (see, e.g., [2l[; [2^; [3^;[i3|)- Outside this 
layer, ideal MHD can be employed to describe the plasma 
motion. The dissipative layer is treated as a surface of 
discontinuity. In order to solve Eq. (|37p boundary condi- 
tions are needed for the variables u and P at this surface 
of discontinuity. In linear theory these conditions are 
described by the explicit connection formulae that deter- 
mine the jumps in the quantities u and P. In order to 
derive the nonlinear counterpart of connection formulae 
we first define the jump of a function, f{x), across the 
dissipative layer as 



[/]= lun{f{x)-f{~x)}.. 



(71) 



where the coordinate system has be translated such that 
Xc = 0. The thickness of this dissipative layer, Sc, is 
determined by the condition that the first and third terms 
in Eq. are of the same order, i.e. 



V^kX 



(72) 



where k is the wavenumber of waves. It is instructive to 
introduce a new, dimensionless variable, a = x/5c, in the 
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dissipative layer. Let Xq be the characteristic width of the 
overlap regions of the dissipative layer (where both the 
linear ideal MHD equations and the nonlinear dissipative 
MHD equations are valid). One of the main reasons we 
have introduced the variable a is the property that a — 
in the dissipative layer, while |a;| — > xo corresponds 
to I (7 1 — > oo. This provides us with the second definition 
of the jump in the function f{x) across the dissipative 
layer, 



[/]= hin {/(a) 

(T—^-tOO 



(73) 



The first connection formula can be obtained in a 
straightforward way by taking into account that the vari- 
able P*^^^ does not change across the dissipative layer, so 
there cannot be any jump in the total pressure. 



[P] = 0. 



(74) 



This connection formula is the same as obtained previ- 
ously by linear and nonlinear theories. 

In order to derive the second connection formula we 
use the approximate relations u w eu^^\ w e^^'^v'^^\ 

P w eP^^^ and introduce the new dimensionless variable, 
g, defined as 



kVSr COS 



a (1) 



(75) 



In the new variable, Eqs. (|56p and (|66p are rewritten as 

V dq 



du 



k cos^ a do ' 



(76) 



dqdq ^ kV^ dP 

da 86 pQv\ |A| dO' ^ ' 



where 



A = R 



2<JA| [(t, + 1)^;2^ + 3c|J 



kV^ 
kV^^ 



(78) 
(79) 



We have used a slightly different Reynolds number, R, 
to the one used in previous sections, here it is defined as 



It is easy to show that the estimations 

q = ff{e^/^kknhR-^), 

A^ff{R^k~Hr^\), 



(80) 



are valid. The ratios of the nonlinear to dissipative term, 
dispersive to dissipative term, and dispersive to nonlinear 
term in Eq. ([77|. N, Dd and £>„ respectively, are 



N = ^(e^^^Rkk. 



(81) 



(82) 



(83) 



The parameters N and Dd can be considered as nonlin- 
earity and dispersive parameters, respectively. Nonlin- 
earity (dispersion) is important if > 1 [Dd !)• When 
N [Dd ^ 1) the nonlinear (dispersive) term in Eq. 
(I77p can be neglected. Dispersion dominates nonlinear- 
ity if Dn > 1. In the opposite case, I?„ < 1, nonlinearity 
dominates dispersion. When klinh = ^(1) (in fa-ct our 
analysis is valid when e^/^ <C kUnh -C e^^^^), the total 
Reynolds number, R, used in this section is of the same 
order of magnitude as that used in the previous sections, 
and the criterion of nonlinearity coincides with that ob- 
tained in section II from the qualitative analysis. With 
the above scalings in mind, it is obvious that the physical 
background of this paper is applicable for relatively short 
inhomogeneity scales. 

Following linear studies, the outer solution reveals that 
v\\ = ff{x~'^) as X 0. Thus, to match the outer and 
inner solutions in the overlap regions the asymptotic re- 
lation q = ff{a~^) as |cr| oo must be valid. It then 
directly follows from Eq. (I77p that 



kV^Pcje) 



(84) 

for \a\ oo. From Eqs. (|76|) and (|84p . we obtain that 

.4 ^„o2. 



Vc% cos a dPr 1 



(85) 



for a ±oo. The functions u+ and w_ are related by 



V 



k cos^ 



dq 
80 



da. (86) 



Equation (|86p uses the symbol of Cauchy principal part, 
because the integral is divergent at infinity. So, in 
accordance with Eq. (j73p . we obtain the implicit jump 
condition for the normal component of velocity 



V 



k cos^ a 



(87) 



Equation (|87p is the nonlinear analog for the implicit 
connection formula for the normal component of velocity. 
The main difference between the linear and the nonlinear 
connection formula is that while in the linear version the 
jumps are expressed explicitly in terms of P{0) and equi- 
librium quantities, in the nonlinear version the jump in 
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the normal component of velocity, u, is expressed implic- 
itly in terms of an unknown quantity q. The connection 
formula ([87| is identical to that found by Ballai et al. 
[25I and Rudcrman et al. 27 1 in the limit of weak non- 



linearity for non-dispersive plasmas. To find solutions in 
the dissipative layer we have to use Eqs. (f37|) and ([77|) 
simultaneously. The boundary conditions for the outer 
solution are provided by Eqs. ((7^1) and ((57)) . 

The governing equation for g, Eq. ([77]) . differs to that 
derived by Ballai et al. [25| only in the appearance of 



the last term in the left-hand side which is related to the 
consideration of the Hall term in the generalized Ohm's 
law. 



V. CONCLUSIONS 

In the present paper we have further developed the 
nonlinear theory of resonant slow MHD waves in the dis- 
sipative layer in one-dimensional planar geometry in plas- 
mas with strongly anisotropic viscosity and thermal con- 
ductivity by considering dispersive effects. The plasma 
motion outside the dissipative layer is described by the 
set of linear, ideal MHD equations. This set of equa- 
tions can be reduced to Eq. ([57)1 for the component of 
the velocity in the direction of the inhomogeneity, u, and 
the perturbation of total pressure, P. The wave motion 
in the dissipative layer is governed by Eq. l|77|) for the 
quantity g, which is the dimensionless component of the 
velocity parallel to the equilibrium magnetic field, de- 
fined by Eq. (j75p . The dissipative layer is considered 
as a surface of discontinuity when solving Eq. ([57]) to 
describe the wave motion outside the dissipative layer. 
The jumps across the dissipative layer are given by Eqs. 
([7^ and ([57]) , thus providing the boundary conditions at 
the surface of discontinuity. In stark contrast to linear 
theory, the jump in u is not solvable analytically, as it 
given in terms of a infinite integral of g - which in turn 
is determined by Eq. (j77p . Since this equation has not 
been solved analytically we must solve Eqs. ([57]) and (1771) 
simultaneously when studying resonant slow waves that 
are nonlinear in the dissipative layer. 

The reader should note that in the upper solar corona 
the plasma (3 (ratio of kinetic to magnetic pressure) is 
very small, hence the significance of slow resonance is 
dramatically reduce (in the limit /3 = 0, slow waves cease 
to exist). The best applicability for the present paper is 
in the regions of the chromosphere and lower solar corona. 

It is interesting to note that the dispersion at the slow 
resonance appears in the form of a nonlinear term [non- 
linear dispersion). It is expected that if the Hall term is 
included in the Alfven resonance, this will appear as a 
linear term. The governing equations and the jump con- 
ditions will be used later in studying the absorption of 
an external driver in the limit of weak and strong nonlin- 
earity. It is intended that the authors will investigate the 
possibility of describing the propagation of solitary and 
shock waves in the slow dissipative layer. It remains to 



be investigated how the inclusion of the dispersive term 
will influence the conclusions drawn when resonant ab- 
sorption was used to explain the damping of waves and 
oscillations in the coronal seismology framework. 



Earlier studies (see, e.g. [16|; [4l|), show that nonlin- 
earity decreases the effective absorption of waves in the 
limit of weak nonlinearity. The effect of dispersion on 
the absorption of waves in slow dissipative layers will be 
addressed in the near future. 
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APPENDIX: THE DERIVATION OF THE HALL 
TERM IN THE INDUCTION EQUATION FOR 
SLOW RESONANT WAVES IN THE 
DISSIPATIVE LAYER 

In this appendix we will derive the components of the 
Hall term in the induction equations and study the condi- 
tions under which this extra effect is important. The par- 
allel component of the magnetic field perturbation domi- 
nates the other components in the slow dissipative layer. 
The Hall term contains the first derivative of this paral- 
lel component of the magnetic field perturbation with re- 
spect to z. The first term of Braginskii's viscosity tensor 
contains the second derivative of the parallel component 
of the magnetic field perturbation with respect to z. As 
a result the Hall term can be of the same order or larger 
than the dissipative term. 

The generalized Ohm's law including the Hall term can 
be written as (see, e.g., |29| ) 



E = -V X B + -j + — j X B, 



(Al) 



where E is the electric field, j the density of the electrical 
current, Uf. the electron number density, e the electron 
charge and a the electrical conductivity. The density of 
electrical current and magnetic induction, B, are related 
by Ampere's law 



j = -VxB, 
with the electrical conductivity given by 



a = 



(A2) 



(A3) 
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Here rUe is the electron mass, Tg the electron collision 
time and r„ the neutral collision time. 

For a fully-ionized, collision-dominated, plasma Eq. 
i3\\ reduces to 



(A4) 



In accordance with Spitzer [42] the electron collision time 
is given by 



2^3/2 

Te = 2.66 X 10^ — ^s, 
n„ In A 



(A5) 



where T is the temperature and In A is the Coulomb log- 
arithm (here taken to be 22). From Eq. (|A5|) . changes 
from 9.4 x lO^^s in the upper photosphere to 1.4 x 10~^s 
in the solar corona. On the other hand, uje changes from 
1.8 X 10^°s~^ in the upper photosphere to 1.8 x 10®s~^ in 
the solar corona. As a consequence the Hall parameter, 
ujeTe, changes from 1.69 x 10^ in the upper photosphere 
to 2.52 x 10^ in the solar corona. Since oJeTe ^ 1, the 
Hall term cannot be neglected in the upper photosphere 
nor the solar corona. 

In order to estimate the relative importance of the Hall 
term and viscous term in the dissipative layer we must 
employ a more sophisticated analysis similar to the anal- 
ysis presented by Ruderman et al. (27j . The generalized 
induction equation (neglecting finite electrical resistivity) 
is. 



— = V X (v X B) + — V X I — B X V X B 

at /le \ne 



(A6) 



In what follows we assume that the ionization coefficient 
is constant, so that is proportional to p, and in partic- 
ular n^^Wue = p~^Wp. Equations ([5^ and (|5D|) provide 
the following estimations in the dissipative layer: 



u 



P = 



(A7) 



where e still denotes the dimensionless amplitude of os- 
cillations far away from the dissipative layer. 

The thickness of the dissipative layer divided by 
the characteristic scale of inhomogeneity is Sc/hnh = 
£?(e^/^). This gives rise to 



dh _ 
' dx 



dh 

'' dz 



(A8) 



where h denotes any of the quantities w, p, v\\ or B\\ . 
Since the first term in the expansion of is independent 
of x, it follows that 



dBx_ 
dx 

d^B± 
dx^ 



dB, 



J, 



= ^(e-i/2ij_j_), 



(A9) 



(the same applies to the variable v±). 

We now need to calculate the components of the vec- 
tors of the Braginskii's viscosity tensor and the Hall term 
from Eq. (|A6[1 normal to the magnetic surfaces (the x- 
dircction) and in the magnetic surfaces parallel and per- 
pendicular to the equilibrium magnetic field lines. We 
use Eqs. (jA8P and (|A9P in order to estimate all the terms 
and we only retain the largest. The components of the 
Braginskii tensor acting in the normal and perpendicu- 
lar directions relative to the equilibrium magnetic field 
are the second and third ones, describing shear viscosity 
(even though in the paper we only consider Eq. ^ we 
need the second and third terms of Braginskii's viscos- 
ity tensor to complete our scalings). Since they are of 
the same order, for the purpose of our estimations it is 
enough to consider rji only. Braginskii's viscosity tensor 
simplifies to, 



??i(V-5i), = r;i — 



77i(V • 5i)_L = Vi 



dx'^ 



^(V • S'o)j| = ?7ocosa ^ 



2 cos a' 



dz^ dxdz 



(AlO) 
(All) 

(A12) 



It should be stated that 77Q ^ ri^ and %(V • 5*1)2: = 
77o(V • S'i)j^ — 0. The components of the Hall term in the 
induction equation reduce to 



H. = 
Hi_ = 
Hi 



Bn COS a sin a d^B 



Bn COS a d'^Bn 

— + . . 

fieue dxdz 

Bq sin a dB\\ dp 
po^eue dz dx 



(A13) 
(A14) 
(A15) 



With the aid of Eqs. (|X7)) . (|X8|) and (|X9|) we obtain 
the ratios 



%(V-5i)x 
'7o(V • -5*0)11 PoVo' 



Vi 

X 



(A16) 
(A17) 
(A19) 



Where x = rjujeTe is the coefficient of Hall conduction and 
rj = 1/ufjL is the magnetic diffusivity. Strictly speaking, 
even the diffusivity is anisotropic in the solar corona, but 
the parallel and perpendicular components only differ by 
a factor of 2. It has been noted that magnetic diffusivity 
is much much smaller that the compressional viscosity 
in the solar corona. However, in the coefficient of Hall 
conduction, x = V^eTei we observe that the magnetic 
diffusivity is multiplied by the product LOgTe, which is 
very large in the solar corona (10^ — 10^). Moreover, if 
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we look at Eq. (|A19[) we see that the coefficient of Hall 
conduction is divided by the density, which is very small 
under solar coronal conditions. Therefore, the parallel 
component of the Hall term in the induction equation 
becomes very important in the slow dissipative layer. 

The Hall terms in the normal and perpendicular di- 
rection relative to the background magnetic field are in- 
cluded for completeness in the paper, but they do not 
play a role in the governing equation (i.e. can be left 
out completely and will not alter the result shown) . This 



is attributed to the fact that the dominant dynamics of 
resonant slow waves is in the parallel direction relative 
to the ambient magnetic field. 

In summary, the Hall term in the parallel direction rel- 
ative to the ambient magnetic field, , must be included 
in the slow dissipative layer when oJeTe 3> 1 because it is 
the same order of magnitude (or larger) than the com- 
pressional viscous term. 



J. Tataronis and M. Grossmann, Z. Phys 261, 203 (1973). 
W. Grossmann and J. Tataronis, Z. Phys 261, 217 
(1973). 

L. Chen and A. Hasegawa, Phys. Fluids 17, 1399 (1974). 

A. Hasegawa and L. Chen, Phys. Fluids 19, 1924 (1974). 

J. A. lonson, Astrophys. J. 226, 650 (1978). 

L. Ofman and J. M. Davila, J. Geophys. Res. 100, 23427 

(1995). 

A. J. C Belien, P. C. H. Martens and R. Keppens, Astro- 
phys. 526(1), 478 (1999). 

M. S. Ruderman and B. Roberts, Astrophys. J. 577, 475 
(2002). 

M. Goossens, J. Andries and M. J. Aschwanden, A&A 
394, L39 (2002). 

J. Terradas, I. Arregui, R. Oliver, J. L. Ballester, J. An- 
dries and M. Goossens, Astrophys. J 679, 1611 (2008). 
M. S. Ruderman, D. Berghmans, M. Goossens and S. 
Poedts, A&A 320, 305 (1997). 

M. S. Ruderman, M. Goossens, J. L. Ballester and R. 
Oliver, A&A 328, 361 (1997). 

W. J. Tirry, D. Berghmans and M. Goossens, A&A 322, 
329 (1997). 

J. M. Davila, Astrophys. J. 317, 514 (1987). 

R. Erdelyi, Sol. Phys. 171, 49 (1997). 

M. S. Ruderman, M. Goossens and J. V. HoUweg, Phys. 

Plasmas 4, 91 (1997). 

S. Poedts, M. Goossens and W. Kerner, Solar Phys. 123, 
83 (1989). 

S. Poedts, M. Goossens and W. Kerner, Astrophys. J. 
360, 279 (1990a). 

S. Poedts, W. Kerner and M. Goossens, Comput. Phys. 
Commun. 59, 75 (1990b). 

S. Poedts, M. Goossens and W. Kerner, Comput. Phys. 
Commun. 59, 95 (1990c). 

T. Sakurai, M. Goossens and J. V. HoUweg, Solar Phys. 



133, 227 (1991). 
[22] J. V. HoUweg, J. Geophys. Res. 90, 7260 (1985). 
[23] J. V. HoUweg and G. Yang, J. Geophys. 93, 5423 (1988). 
[24] M. S. Ruderman and M. Goossens, Astrophys. J. 471, 

1015 (1996). 

[25] I. Ballai, M. S. Ruderman and R. Erdelyi, Phys. Plasmas 

5(1), 252 (1998). 
[26] M. S. Ruderman, J. Plasma Phys. 63, 43 (2000). 
[27] M. S. Ruderman, J. V. HoUweg and M. Goossens, Phys. 

Plasmas 4, 75 (1997). 
[28] S. I. Braginskii, Rev. Plasmas Phys. 1, 205 (1965). 
[29] E. R. Priest, Solar Magnetohydrodynamics (D. Reidel, 

Dordrecht, 1982). 
[30] R. Erdelyi and M. Goossens, A&A 294, 575 (1995). 
[31] V. B. Baranov and M. S. Ruderman, Fluid Dyns. 9, 

(1974). 

[32] M. S. Ruderman, Fluid Dyns. 11, 762 (1976). 
[33] M. S. Ruderman, Fluid Dyns. 22, 229 (1987). 
[34] M. S. Ruderman, J. Plasma Phys. 67, 271 (2002). 
[35] I. Ballai, J. C. Thelen and B. Roberts, A&A 404, 
(2003). 

[36] R. Miteva, I. Zhelyazkov and R. Erdelyi, New J. Phys. 

6, 14 (2004). 
[37] J. D. Huba, Phys. Plasmas 2, 2504 (1995). 
[38] A. H. Nayfeh, Introduction to Perturbation Techniques 

(Wiley-Interscience, New York, 1981). 
[39] M. Goossens, M. S. Ruderman and J. V. HoUweg, Solar 

Phys. 157, 75 (1995). 
[40] J. V. HoUweg, Astrophys. J. 335, 1005 (1988). 
[41] I. Ballai, R. Erdelyi and M. S. Ruderman, Phys. Plasmas 

5(6), 2264 (1998). 
[42] L. Spitzer, Physics of Fully Ionized Gases (Interscience, 

New York, 1962). 



421 



701 



